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GAS-RICH MERGERS IN LCDM: DISK SURVIVABILITY AND THE BARYONIC ASSEMBLY OF GALAXIES 
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ABSTRACT 

We use A'^-body simulations and observationally-normalized relations between dark matter halo mass, stellar mass, 
and cold gas mass to derive robust expectations about the baryonic content of major mergers out to redshift z ~ 2. 
First, we find that the majority of major mergers [m/M > 0.3) experienced by Milky Way size dark matter halos 
should have been gas-rich, and that gas-rich mergers are increasingly common at high redshift. Though the frequency 
of major mergers into galaxy halos in our simulations greatly exceeds the observed early-type galaxy fraction, the 
frequency of gas-poor major mergers is consistent with the observed fraction of bulge-dominated galaxies across the 
halo mass range Mdm ~ 10^^ — 10^"^ M©. These results lend support to the conjecture that mergers with high baryonic 
gas fractions play an important role in building and/or preserving disk galaxies in the universe. Secondly, we find that 
there is a transition mass below which a galaxy's past major mergers were primarily gas-rich and above which they were 
gas poor. The associated stellar mass scale corresponds closely to that marking the observed bimodal division between 
blue, star-forming, disk-dominated systems and red, bulge-dominated systems with old populations. Finally, we find 
that the overall fraction of a galaxy's cold baryons deposited directly via major mergers is significant. Approximately 
~ 20 — 30% of the cold baryonic material in Mgtar ^ lO^^'^Af© (A^dm ^ lO^^M©) galaxies is accreted as cold gas or 
stars via major mergers since z = 2, with most of this accretion in the form of cold gas. For more massive galaxies 
with Mstar ^ IO^^Mq (Mdm ^ IO^^Mq) the fraction of baryons amassed in mergers since z = 2 is even higher, ~ 40%, 
but most of these accreted baryons are delivered directly in the form of stars. This baryonic mass deposition is almost 
unavoidable, and provides a limit on the fraction of a galaxy's cold baryons that can originate in cold flows or from 
hot halo cooling. 

Subject headings: cosmology: theory — dark matter — galaxies: formation — galaxies: halos — 
methods: A^-body simulations 



1. INTRODUCTION 



In the cold dark matter (CDM) model of structure 
formation, major galaxy mergers are believed to play 
an impor tant role in determining a galaxy's morphol- 
ogy (e .g. iToomre &: Toqmrd 11972: iBarnes fc Hernguistl 



[1991 [Robertson et all l2006al lbl: iBurkert et all l2008f ). 
well as triggering star formation and AGN 



as 



ac- 



tivity fe.g. | Mihos fc Hcrnaiusti 11996: iHeckman et al.l 

'laOOSD, while mi- 



tivity I e.g. iivimos &: ncrnamsti liijiJ 
[19861 ISpringel et all 120051: ICox et all 



nor mergers may help explain the origin of thick disks 
and extended diffuse lig ht cor nponents around galax- 
ies (IBarnes fc Hernouistl 119961: iK azantz idis et all 120081: 
lYounger et al. 2007nPurcell et all 
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l2008HViUalobos fc Hclmi 2008: ' Kazantzidis et al.ll2009n . 
More than simply triggering star formation in existing 
gas and altering existing galaxy morphologies, mergers 
deliver new stars and additional fuel for star formation, 
and thereby contribute to baryonic acquisition of galax- 
ies over their histories. That mergers contribute sig- 
nificantly to many aspects of galaxy formation is now 
fairly well accepted, however there are lingering concerns 
that mergers are too common in CDM to explain the 
prominence o f thin disk-dominated galaxies in the local 
universe ("e.g. lToth fc Ostri keil[r992|: IWalker et al.lll996l: 

nL2 " ' ' ' 



These sim- 
mass scale 



Stewart et aLl 72008: Purc cU et al.l l2009t iBullock et all 
20091 and references therein). Here we explore the bary- 



onic content of these predicted mergers and the potential 
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ramifications of gas-rich and gas-poor mergers on galac- 
tic morphological evolution. 

The baryonic delivery of material into galaxies via ma- 
jor mergers touches on a broader question in galaxy 
formation: how do galaxies get their baryons? In re- 
cent years, studies motivated by hydrodynamic simula- 
tions have placed a growing emphasis on the importance 
of smooth gas accretion via "cold flows." These cold 
flows constitute streams of cold gas flowing along fila- 
mentary structures (particularly at high redshift) with 
sufficiently high densities to penetrate into a halo's cen- 
tral region w ithout heating the gas to the virial tem- 
perature (e.g.lBirnboim fc I3 ekol 2003: Keres et al.ll2005l: 
Dekel fc Birnboimlj200^ nKcres et al. 200a; ,Dekel et al l 
2009HBrooks et al.|[2009l: [Agcrtz et al. 2003). 
ulations demonstrate a characteristic halo 
(~ IO^^Mq) below which cold streams are the dominant 
mode of gas accretion, and above which gas cooling di- 
rectly from shock-heated (hot mode) material dominates. 

Though there has yet to be any observational evidence 
that cold fiows actually occur in nature, the possibility is 
well motivated by theory and is suggestive of a number 
of interesting scenarios for galaxy assembly. One partic- 
ularly interesting idea is that flows of cold gas are vital 
to the formation o f disk galaxies at high redshift z > 1 
(IDekel et al.ll2009l) . StiU, even if disks were built at high 
redshift via streams of cold gas, we can return to the issue 
of disk survival raised above. How do obs erved popula- 
tions of disk galaxies at high redshift (e.g. I Wright et al.l 
at z ~ 1.6) survive subsequent mergers and remain 
disk-dominated hy z = 0? 

In a previous paper ([Stewart et al.l I2008D we studied 
the merger histories of Milky Way-size dark matter ha- 
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los within a cosmological iV-body simulation and found 
that approximately 70% should have accreted an ob- 
ject with more than twice the mass of the Milky Way 
disk (m > IO^^Mq) in the last 10 Gyr. In order to 
achieve the ~ 70% disk-dominated fraction that has been 



observed in Milky Way-sized halos ( 


Weinmann et al. 
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;lPark et al.ll2007HChoi et al.ll2007l: 


Weinmann et al. 


2008 


), mergers involving > 1/3 mass-ratio events must 



not always destroy disks. Adding to the associated 
concern, Purcell, Kazantzidis, and Bullock (2008b) per- 
formed focused numerical experiments to study the im- 
pact of m = IO^^Mq encounters onto fully-formed Milky- 
Way type thin stellar disks and concluded that thin 
(~ 400 pc) dissipationless stellar disks do not survive 
these (presumably common) encounters. 

One possible solution appeals to the role of cold gas 
in mergers. Focused merger simulations in the past few 
years have begun to suggest that sufficiently gas-rich 
mergers may help build angular momentum in the cen- 
tral galaxy, while feedback physics prevents the gas from 
forming stars too quickly, r esulti n g in a disk-dominated 

[200l iSoringel fc HcrnauistJ 



merger remnant (IBarne! 



20051 iRobertson et al.[ l2006at iBrook et al 2007b|: 
Hopkins et all l2009al : IRobertson fc Bullock! i200i ). 

Encouragingly, cosmological simulations that have 
been successful in reproducing disk galaxies have 
also shown that gas-rich mergers have pl ayed an im- 
portant role in the d i sk's c r eation (e.g. iBrook et al.l 
120041 : iGovernato et all 120071 . 12009( 1 . and there have 
been recent observations of late-type galaxies that 
may be in the process of re forming after a recen t 

gas-rich major mer ger (e.g. iHammer et al.l l2009b[ ). 

IRobertson fc BuUockl (|2008[ ) also showed that the disk- 
like merger remnants from gas-rich mergers are similar 
to the kine matically hot disks that have b een observed 
at z ~ 2 (iForster Schreiber et all 120061 : iGenzel et all 
l2006t IShapiro et al.ll2008fr 

Our goal in this work is to provide an empirically- 
motivated accounting for the expected gas and stellar 
content of mergers by relying on robustly determined 
dark matter halo merger rates. We aim first to determine 
whether gas-rich mergers are common enough to signifi- 
cantly alleviate the disk formation problem (a necessary 
but not sufficient condition in evaluating this scenario). 
We also aim to investigate the overall importance that 
mergers play in the acquisition of a galaxy's cold baryons. 

While the merger histories of the dark matter halos 
are predicted accurately in dissipationless A^-body simu- 
lations, the baryonic content of these mergers arc much 
more difficult to predict from first principles. Indeed, 
an accurate ab initio accounting of the baryonic con- 
tent of dark matter halos would require an overarching 
theory that solved all of the major problems in galaxy 
formation, including star formation, feedback, and the 
complicated interplay between mergers and galaxy as- 
sembly. In this work we chose to avoid the issues of 
galaxy formation physics entirely. Instead we adopt a 
semi-empirical approach that forces our model to match 
observations at various redshifts. First, we adopt the 
technique of monotonic abundance matching to assign 
stel lar masses to dark matte r halos (specifically follow- 
ing |Conrox&!^^chsli3 HOOl) . Second, we use observa- 
tional relations between stellar mass and gas mass (e.g. 
lMcGaugh.2005l : lErb et al.„200Gl to assign gas masses to 



our halos. We then combine these relat ions with the 
- /V-bod y hal o merger histo r ies de scribed in lStewart et al.l 
(|2008f) and iStewart et al] (|2009f) in order to determine 
the baryonic properties of mergers back to redshift z ~ 2. 
In §2 we discuss the details of our method. Wc present 
our primary results and discuss the implications for disk 
survival and the baryonic assembly of galaxies via merg- 
ers in §3. We summarize our main conclusions in §4. 

2. METHOD 
2.1. The Simulation 

Our simulation contains 512'^ particles, each with mass 
nip = 3.16 X IO^Mq, evolved within a comoving cu- 
bic volume of 80h~^ Mpc on a side usin g the Adaptive 
Refinement Tree (ART) A^-body code ( Kravtsov et al.l 
[19971 12OOI . We assume LCDM cosmological parame- 
ters: = 1 - f^A = 0.3, h = 0.7, and erg = 0.9. The 
simulation root computational grid consists of 512'^ cells, 
which are adaptively refined to a maximum of eight lev- 
els, resulting in a peak spatial resolution of 1.2h~^ kpc 
(comoving). We give only a brief overview of the sim- 
ulation here, as it has been reviewed mo re extensively 
in pre vious works. We refer the reader to IStewart et al.l 
(|2008f ). and reference therein, for a more complete dis- 
cussion. 

Field dark matter halos and subhalos arc identified us- 
ing a variant of th e bound density maxima algorithm 
(jKlvpin et al.lll999l ) . A suhhalo is defined as a dark mat- 
ter halo whose center is positioned within the virial ra- 
dius of a more massive halo, whereas a field halo does 
not. The virial radius is defined as the radius of a col- 
lapsed self gravitating dark matter halo within which the 
average density is Avir times the mean density of the 
universe. Under comparison to constructed mass func- 
tions, we have determined that our halo catalogs are com- 
plete to a minimum mass of IO^'^A/q, and our sample 
includes ^ 17, 000 field halos at z = in the mass range 
M = 1011- 2-13. 2j^^^^ -yyg usc th e Same merger trees de- 
scribed in IStewart et al.l ( 120081 ). constructed using the 



1 (120081). constructed using 
Wechsler et all (12^^ . 12001 . 



techniques described in 

Wc present our results primarily in terms of the dark 
matter mass ratio between the two halos that are un- 
dergoing a merger, (m/M)DM, where we always define 
mDM as the mass of the smaller dark matter halo (which 
we will sometimes refer to as the satellite halo) just 
prior to entering the virial radius of the larger one, 
and Mdm is the mass of the larger dark matter halo 
(also referred to as the host halo) at this infall epoch. 
Thus, Mdm does not incorporate the mass moM, and 
(m/M)DM has a maximum value of 1.0. However, we 
also present results in terms of the stellar mass ratio of 
the central galaxies within merging halos, or the mass 
ratio between the total baryonic mass of these central 
galaxies (stellar mass plus gas mass). We refer to the 
dark matter, stellar, and galaxy (baryonic) mass ratios 
as (m/M)DM, (TO/M)star, (n^/M)gai, respectively. Inde- 
pendent of the mass ratio definitions above, we always 
refer to a merger ratio of m/M > 0.3 as a major merger. 
For the sake o f comparison to our past work, we em- 
phasize that in lStewart et al.l (|2008f ). we considered two 
definitions of merger ratio. The first (written as to/Mq in 
that paper) referred to the ratio of the satellite mass at 
infall m to the final dark matter halo mass Mq at z = 0. 
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Fig. 1. — Two step method for assigning ba ryons to dark matter halo s. Left: Stellar mass Mstar versus dark matter halo mass Mdmi 
for 2 = 0, 1,2, based on abundance matching l lConrov Ik We chsler 20091) . Middle(Right): Power law fits to Mgas/Mgtar as a function of 
stellar mass at z = 0(2), wit h the correspond ing gas fraction /gas = A-fgas / ( Mpas + Afatar ) sho wn on the right axis. The symbols represent 
observational estimates from lMcGaughll2005l (entire sample: blue triangles) . Ifvannappanir2004l (av erage of blue galaxy sa mple, gold squares; 
average of re d galaxy s ample, gold X's). lEaldrv et al.ll200^ (average: green short-dashed line). 'Somcrvi lle et al.ll2008l (average: red dot- 
dashed line). IWei et al.^ 2009 (best fit to total sample: magenta dot-dot-dot-dashed line) and Erb et al. 20061 (entire sample: red triangles). 
The solid (black) line in each panel is the best-fit relation (Equation [TJ, while the long-dashed and dotted (black) lines represent the la 
and 2(T scatter, respectively. 



This is not the ratio we are using here. The mass ra- 
tio definition we adopt here is more standard and refers 
always to the mass ratio at the redshift z of accretion 
[m/Mz in the notation of Stewart et al. 2008). 

In order to provide robust results, we define a merger 
to occur once the smaller halo crosses within the virial 
radius of the larger halo and becomes a subhalo. as 
the subsequent orbital evolution of each subhalo will 
depend on the baryonic distribution within both ha- 
los. We emphasize that for the major mergers we 
consider in this paper, the dynamical friction decay 
timescales are expected to be short (comparable to the 
halo dynamical timescale) for t ypical orbital parame- 
ters (jBovlan-Kolchin et al.ll2008[ ). Since only ~ 5% of 
IO^^Mq halos have experienced a major merger in the 
past halo dynamical time at z = (see Figure 2]), most 
of these major mergers into the virial radius have had 
adequate time to impact the central galaxy and do not 
survive as distinct substructure by z = 0. For a more in- 
depth comparison between merger rates into the virial 
radius and the rate at which accreted satellites are "de- 
stroyed" in this simulation (i.e. once the y lose 90% of 
their infall mass) we refer the reader to iStewart et al.l 

2.2. Assigning Stars 

One particularly simple, yet surprisingly success- 
ful approach for assigning galaxies to dark mat- 
ter halos is to assume a monotonic mapping be- 
tween dar k matter halo mass MpM and galaxy lumi- 
nositv L (iKravtsov et~al] 120041: iTasitsiomi et all 120041: 



Vale fc Ostrikeri 120041: iConrov et all 120061: iBerrier et all 
20061: IPurcell" et al. [ I2007t IMarm et all 120081)" Using this 
technique, provided we know ng{> L) (the cumulative 
number density of galaxies brighter than L) we may de- 
termine the associated dark matter halo population by 
finding the halo mass above which the number density 
of halos (including subhalos) matches that of the galaxy 
population, nh{> Mdm) = "■g(> L). 

We use a similar approach, and instead assume a 
monotonic relationship between halo mass and stellar 



m a.ss Mstar- Specifica l ly, we adopt the relation found 
bv lConrov fc Wechsleri (|2009t ) (hereafter CW08; interpo- 
lated from their data as shown in their Figure 2). Figure 
[T] (left panel) shows the resulting relation between stellar 
mass and dark matter halo mass for z = 0, 1,2 (upper 
blue, middle black, lower red curves respectively) where 
Mstar is the stellar mass of the central galaxy residing 
within a dark matter halo of mass Mdm- We ignore scat- 
ter in this relationship for the results that follow, but we 
find that including a Gaussian scatter of 0.1 dex has no 
substantive effect our results. 

Of course, a simple relation of this kind cannot be cor- 
rect in detail, however, in an average sense, it provides 
a good characterization of the relationship between halo 
mass and galaxy stellar mass that must hold in order for 
LCDM to reproduce the observed universe. Moreover, 
by adopting it wc insure that our model self-consistently 
reproduces the observed stellar mass function of galaxies 
out to z ~ 2. We cannot use this method to explore 
merger rates as a function of stellar mass beyond z ^ 2 
because the stellar mass function is poo rly constrained at 
higher redshifts. Wc refer the reader to lMarchesini et all 
(j2009f l for a detailed investigation of random and sys- 
tematic uncertainties in computing the stellar mass func- 
tion at 1.3 < z < 4.0 (e.g. the impact of different SED- 
modeling assumptions, cosmic variance, and photometric 
redshift errors). 

2.3. Assigning Gas 

In order to reasonably assign gas to the central galaxies 
within our halos, we quantify observationally-inferred re- 
lations between the ratio of cold gas mass to stellar mass 
(Mgas/Mstar) as a function of st ellar mass using the em- 
pirical results of lMcGaughll2005l (bl ue triangles, for disk- 
dominated galaxies at z = 0) and lErb et all 120061 (red 
triangles, for UV-selected galaxies at z 2), as shown in 
the middle and right panels of Figure [T] Though both of 
these samples are biased with respect to blue (gas-rich) 
galaxies, we argue below that by adopting these relations 
we are not strongly biasing our overall results. 

As shown by the black solid lines in the right panels of 
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Figure [H the z ~ and z ^ 2 cold gas fraction data can 
be characterized by a relatively simple function of stellar 
mass and redshift: 



^ - 0.04 



4.5 X 10"Af(7 



-a(z) 



(1) 



where the gas fraction relation evolves and steepens 
with redshift as a[z) ~ 0.59(1 + z)^'^^^ . Assuming a 
Gaussian scatter about the best-fit lines (independent 
of mass), we find that the scatter evolves with redshift 
as logiQ [cr(2;)] = 0.34 — O.lOlog^g (1 -I- z), such that the 
correlation between the cold gas fraction and stellar mass 
is tighter at z ^ 2 than it is at z = 0. The black long- 
dashed and dotted lines in the right panels of Figure 1 
demonstrate the la and 2cr scatter, respectively. In order 
to assign gas content to our halos as a function of stellar 
mass, we draw randomly from a Gaussian distribution 
with the average value and standard deviation given by 
the above analytic characterizations of Mstar(-Afgas)- 

For comparison, the middle panel of Figure [l] also 
shows several additional observational estimates for the 
gas-fraction relation as a function of stellar mass. The 
gold squares and gold X' s prese nt the average relation 
measured by iKannappanI (|2004f ) for blue galaxies and 
red galaxies, respectively. The green short-dashed line 
shows the average (statistical) relation derived using a 
combination of published galaxy stellar mass functions 
and the ob served stellar mass-metallicity relation by 
IBaldrv eFal.. (i2008il , who used robust chemical evolution 
arguments to derive implied gas fractions as a function of 
stellar mass. Finally, th e avera ge results from direct mea- 
surements bv lWei et ail (|2009D are shown by the magenta 
dot-dot-dot-dashed line. For the sake of comparison, we 
sh ow the predicted relati on from the semi-analytic model 
of lSomerville et aLll2008l (red dot-dashed line) . While we 
do not utilize these additional data sets in constructing 
our fitting function, they conform well to our average re- 
lation and certainly lie within the Icr scatter of our fit 
to the McGaugh ( 2005) data (with the ex ception of low- 
mass red galaxies from iKannapparJ 120041 : see discussion 
below). We also not e that the cold gas fractions de- 
rived bv lWright et al.l (|2009[ ) for six galaxies at z ~ 1.6 
are also consistent with the evolution in our fit, with 
every galaxy in their sample falling within our 2a scat- 
ter. Their sample does have a slightly higher average 
gas fractions at fixed A/gtar than our adopted relation, 
but the discrepancy is not significant given the small- 
number statistics. We have also compared our f i t to th e 
34 galaxies at z ~ 0.6 studied in lHammer et al.l ()2009af ). 
in which the authors use K-band magnitudes to esti- 
mate t otal stellar mass via the methodology of Bell et al.l 
(l200l . and assume the Schmidt-Kcnnicutt law to derive 
gas masses from star formation rates. Encouragingly, 26 
of their galaxies {r^ 75%) fall within the la contours of 
our best-fit at this redshift, with all 34 of them within 
2a. 

The fact that we have fit our z = relation to disk- 
dominated galaxies introduces a potential worry about 
applying the relation to every galaxy halo in our sim- 
ulation, including ones that presumably host massive 
(spheroidal) galaxies. However, it is unlikely that this 
bias will drastically affect our results, primarily because 
the gas fractions in the adopted relation are only appre- 



ciable (> 0.5), in the smallest galaxies (at z = 0) with 
A'/star S ^O'^'^'^Mq - the stellar mass regime that is known 
to be dominated by disk-dominated galaxies (see, e.g., 
the left panel of Figure 2). For larger galaxies, it is reas- 
suring to no te that the average relation for the red galaxy 
sample from lKannappanI ()2004f) lies within our adopted a 
scatter and is in relatively good agreement with the other 
(disk-selected) observations. It is only for less massive 
galaxies (a regime where blue disk galaxies dominate the 
total population anyway) where the red galaxy sample of 
[Kannap pau ( 20041 ) becomes significantly discrepant from 
our fiducial relation. Finally, even if our fiducial relation 
is biased to be slightly high for massive galaxies, the gas 
fractions are already small enough that we would never 
classify them as "gas-rich" in our discussions below. 

A similar poi nt of concern may be applied to the 
lErb et al.l (|2006f) data at z ~ 2. These galaxies were se- 
lected based on UV luminosity and thus constitute an ac- 
tively star-forming population. However, there is a good 
deal of evidence that UV luminosity is tightly correlated 
with total st ellar mass (or halo mass) at z > 2 (see e.g., 
discussion in lConrov et al.ir2008l . and references therein). 
For example, galaxies with higher U V luminosities at 
z ~ 2 are more strongly clustered (jAdelberger et al.l 
120051 ). suggesting that they reside within more massive 
dark matter halos. In addition, the UV and V-band lumi- 
nosity functions of galaxies at z ~ 3 are in relative agree- 
me nt, producing similar number densities for ^ L* g alax- 
ies (iShaplev et al.ll2001l:ISa"wicki fc Thompsonll2006l : also 

l|l2 



see Table 2, and discussion, in Istewart et all 120091) . As 
such, it is reasonable to consider a galaxy sample selected 
on UV luminosity to contain a fairly representative sam- 
ple of bright galaxies at z ~ 2. 

We also note that the gas estimates we adopt from 
lErb et all (l2006h assume the global Schmidt law of 
iKennicufrB ()1998l ): Esfr oc Sg^t- Both observations and 
recent hydrodynamic simulations have suggested that 
while this relation is tightly correlated for molecular gas, 
it may underestimate the total gas content, especially for 
galaxies whe re the fraction of gas in molecular form is not 
unifo r m (e.g.lWong fc Blitjl2002l : iRobertson fc Kravtsovl 
I200I IGnedin et all 120091). As a conseq uence, the gas 
fraction estimates from " lErb et all (I200I may represent 
a lower limit, such that the evolution of gas fraction with 
redshift may actually be steeper than our adopted rela- 
tion. Insofar as issues of gas accretion and disk survival 
are concerned, our relation may be considered a conser- 
vative lower limit on the estimated gas content of our 
galaxies. 

Finally, there has been some discussion in the liter- 
ature that the stellar initial mass function (IMF) may 
evolve systematically to become more top heavy at 
high redshift in galaxies with extremely low metallicities 
(e.g. 'Lucat ello et al|[2005l : lTumlinson|[2007t Ivan Dokkur3 
[2008; Ko miva et al.ll2008l) . This evolution in the IMF has 
not been corrected for in our adopted mapping between 
halo mass and stellar mass. While we do not expect this 
to significantly affect our results, including this evolution 
of the IMF would decrease our stellar mass estimates at 
fixed halo mass — which would, in turn, increase our es- 
timated gas fractions. As such, so far as issues of gas 
accretion and disk survival are concerned, the results we 
present are a conservative lower limit. 
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With the above qualifications in mind, we now turn 
to the impHcations of this empirically-motivated stellar 
mass and gas mass assignment prescription. 

3. RESULTS AND IMPLICATIONS 
3.1. Galaxy Morphology 

We start by investigating the merger histories of z = 
dark matter halos. The solid black line in the left panel 
of Figure [2] shows the fraction of dark matter halos that 
have experienced at least one major dark matter merger 
with (m/A/)DM > 0.3 since z = 2 as a function of dark 
matter halo mass (lower axis label). Equivalcntly, the 
merger fraction as a function of galaxy stellar mass can 
be seen by focusing on the upper axis label. Compare this 
result to the black squares, which show the early- type 
fraction for SDSS galaxies as a function of central halo 
mass as derived by Weinmann et al. (2006; with "early- 
type" based on galaxy color and specific star formation 
rate) ^. Also compare to the early- type fraction for SDSS 
galaxies as a function of halo mass, where "early-type" is 
defined by the concentration parame ter (C > 3), shown 
as the black crosses, from iWcinman n et al.| |2009) (Wc 
will refer to these two results as W06 and W09, respec- 
tively). Clearly the fraction of halos with major mergers 
greatly exceeds the early-type fraction at low masses. 

Consider now the likely baryonic makeup of these 
mergers. The (blue) dotted line shows the fraction of 
halos that have experienced a gas-rich major merger 
since z = 2 and the (red) dashed line shows the frac- 
tion of halos with at least one gas-poor merger. In our 
fiducial case, we define a merger to be gas-rich if both 
the central galaxy and the infalling satellite galaxy have 
more baryonic mass in the form of gas than in stars: 
fg = Afgas/(Afgas Mstar) > 50%. Similarly, gas-poor 
mergers arc defined such that each of the progenitors has 
fg < 50%. The shading of the red and blue bands corre- 
spond to varying the definition of gas-rich from fg > 30% 
to fg > 70%. 5 

Remarkably, if one makes the simplistic assumption 
that only gas-poor mergers generate early-type galaxies 
(red dashed line) and gas-rich mergers preserve disks, 
then the observed SDSS relation from W06 and W09 

^ Note that in W06, they divide galaxies into three categories 
instead of two; early-type, late-type and intermediate-type. In 
order to compare our simple bimodal model to their findings, wc 
count half of their intermediate-types as early-type, and half as 
late-type. 

Because the morphological impact of a mixed merger (where 
one galaxy is gas-rich and one is gas-poor) is largely unclear, we 
choose to focus on the extreme cases where both are either gas- 
rich or gas poor, and to leave a more detailed exploration of mixed 
mergers for future work. This means that the combined gas-rich 
fractions and gas-poor fractions in Figure [2] need not equal the 
total merger fraction, which includes all major mergers regardless 
of baryonic content. However, because a galaxy's gas fraction is 
a strong function of halo mass at fixed redshift (and because we 
define a strict cutoff between gas-rich and gas-poor based on gas 
fraction) we find that mixed mergers are less frequent than mergers 
between two gas-poor or two gas-rich systems. If the larger galaxy 
in a major merger is gas-rich (by our definition), then the smaller 
galaxy is most likely gas-rich as well. Conversely, if the smaller 
galaxy in a major merger is gas-poor, then the larger galaxy is 
most likely also gas-poor. While there does exist a characteristic 
mass scale for which mixed mergers become a significant portion of 
the overall merger fraction, even at this special mass scale they still 
only constitute about half of all mergers (see Figure |4] discussion 
in §3.3). 



is reproduced fairly well. Specifically, wc find that the 
fraction of halos with a disk-destructive merger increases 
from only - 15% at IO^^Mq to ~ 55% at lO^^Mg, in re- 
markably good agreement with W06 and W09. Not only 
does this agreement provide a possible solution to disk 
survivability within Milky Way-sized halos, but it also 
implies that the gas-poor merger history of a dark mat- 
ter halo may be closely tied to the halo mass-morphology 
relation (across this range in halo mass). Although halo 
merger rates have been shown to depend strongly on en- 
vironment (Fakhouri & Ma 2009), suggesting a possible 
connection to the morphology-density relation of galax- 
ies, we see in Figurc[5]that the overall halo major merger 
rate (solid black line) is not steep enough to account for 
the observed change in morphological fraction with halo 
mass. 

It is also worth mentioning that the implied transi- 
tion between mostly gas-rich and mostly gas-poor merg- 
ers occurs at a characteristic mass Mstar — 5 x IO^^Mq 
(or equivalently Mum — 2 x IQ-^^Mq), which is close 
to the characteristic bimodality scale that separates 
blue, star-forming, disk-dominated systems and red, 
bulge-dominated systems with old populations (typically 
10^°M^see e ' 
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Of course, for detailed treatments of galaxy morphol- 
ogy, this model is too simplistic. The inclusion of gas- 
richness in the efficacy of major mergers to disrupt mor- 
phology seems to greatly relieve the problem of disk sta- 
bility, and the agreement between the observed early- 
type fraction from SDSS and the fraction of halos with 
at least one gas-poor major merger is quite remarkable. 
However, this is only a first step in understanding the 
distribution of galaxy morphologies. In detail, nothing 
we have investigated here can explain the prominence of 
"bulgeless" galaxies, as simulated gas-rich major merg- 
ers lead to galaxies with noticeable bulges and disks that 
are thicker and hotter than the Milky Way. Even cosmo- 
logical simulations th at produce thin disk galaxies (e.g. 
iGovernato eTaI][2009f) require significant smooth gas ac- 
cretion from the hot halo after the most recent gas-rich 
merger in order to form a thin disk. Such intricate details 
are beyond the scope of this paper, as we are primar- 
ily concerned with providing the most robust predictions 
possible, only using A^-body dark matter halo merger 
trees and empirical relations between M-qm, Msta.r, and 
Mgas to determine gas-rich and gas-poor merger statis- 
tics. 

Any model that predicts detailed bulge-to-disk mass 
ratios or estimates the thinness or thickness of the galac- 
tic disk resulting from a merger event must require fur- 
ther assumptions about the detailed morphological ef- 
fects of any given merger event, which are still relatively 
uncertain. We refer the reader to iHopkins et al.l ()2009bl ) 
for a more detailed galaxy formation model that gener- 
ates bulge and disk mass estimates due to merger events 
as a function of merger mass ratio, gas fraction and or- 
bital parameters. Using a semi-empirical assignment of 
gas and stars to dark matter halos similar (but distinct) 
from our own treatment, they find that cosmologically 
motivated merger trees lead to consistent distributions 
of B/T (bulge mass to total galaxy mass) values as a 
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Fig. 2. — Fraction of dark matter halos with at least one major merger since 2 = 2, as a function of host halo mass Mdm (lower axis) 
and stellar mass Mstar(upper axis), for varying definitions of 'major merger.' Left: Major merger defined by the ratio of dark matter halo 
masses, (m/M)p M- The black squares and crosses in this figure show the observed early-type fraction as a function of halo mass from 
IWeinmann et al.l (2006, 2003). Middle: Major merger defined by the ratio of the stellar masses in each central galaxy, {m/Af)star- Right: 
Major merger defined by ratio of the total baryonic mass of the central galaxies, (m/Af)gai. In each panel, the solid (black) line shows 
the total merger fraction. The dashed (red) line shows the merger fraction while only considering mergers between two halos that both 
contain gas-poor central galaxies {fg < 50%). The dotted (blue) line shows only mergers for which both galaxies are gas-rich (/g > 50%). 
The shaded regions surrounding the red and blue lines represent the impact of varying our distinction between gas-rich and gas-poor from 
30% to 70%. The bottom axis shows the same range in halo mass in each panel, while the corresponding stellar (or baryonic) mass of the 
central galaxy is shown on the top axis. Note that the range in y-values in the right panel is larger than the left and middle panels. Error 
bars are Poissonian based on the number of host halos and the total number of mergers, and do not include possible errors in assigning 
stars and gas to halos (though we do account for scatter in the Af star ( Afgas ) relation, see §2.2,2.3). 



function of halo mass and redshift. 

3.2. Alternative definitions for major merger 

The middle and right panels of Figure 2 explore how 
the implied merger fraction trends change when one 
chooses to define major mergers using the stellar-mass 
ratio, (m/M)star > 0.3, and total baryonic galaxy mass 
ratio, (TO/A/)gai > 0.3, respectively, rather than the total 
mass ratio in dark matter. Clearly, the implied trends 
between merger fraction and galaxy halo mass depend 
sensitively on whether dark matter mass ratios, stellar 
mass ratios, or bary onic galaxy ma ss r atios are consid- 
ered (also see lMaller 2008; Stewart 200^. As seen by the 
solid black line in the middle panel, high stellar-mass ra- 
tio events are rare in small galaxy halos and common in 
high mass halos. This follows directly from the fact that 
low-mass halos tend to have a higher stellar-mass to dark 
matter mass ratios (see Figure 1). The trend changes 
dramatically when the full baryonic mass of the galaxy 
is considered in the ratio (right panel). In this case, 
even small galaxy halos are expected to have had com- 
mon mergers with galaxies of a comparable total baryonic 
mass (note that the range of the vertical axis has changed 
in the right-hand panel). It is clear from this compari- 
son alone that most of the major mergers experienced by 
small galaxies must be gas-rich. 

We note that the fraction of systems that have expe- 
rienced at least one gas-rich merger (and consequently, 
the total merger fractions as well) show qualitatively 
different behavior depending on these definitions. Gas- 
rich halo mergers ((m/M)DM > 0.3) are relatively fre- 
quent for Mdm = lO^-^M© systems (40% since z = 2), 
with a smoothly declining merger fraction for increasing 
halo mass (roughly linear in log Mdm), while the gas- 
rich stellar merger fractions decline in a qualitatively 
similar fashion but are universally less common, with 
merger fractions < 20% for Mdm = 10^^ '^ Mq. In con- 



trast, the fraction of halos with at least one major galaxy 
merger ((m/M)gai > 0.3) shows completely different be- 
havior, with extremely high fractions (40 — 90%) and 
non-monotonic evolution with log A/dm (with a maxi- 
mum value at Mdm ^ lO^^-^Afg). 

The behavior of gas-poor merger fractions, on the other 
hand, remains remarkably similar in each case. Regard- 
less of these three merger ratio definitions, the fraction of 
halos which have experienced a gas-poor major merger 
is negligible at small halos masses (Mdm ^ IO^-^'^Mq) 
and increases roughly linearly with log Mdm to a frac- 
tion of 65 - 75% at Mdm = lO^^-^Mg. Because these 
gas-poor merger fractions appear somewhat independent 
of the merger ratio definition used (and they remain con- 
sistent with observed morphological fractions as a func- 
tion of halo mass) we again suggest that a dark mat- 
ter halo's gas-poor merger history may be a particularly 
useful tracer of galaxy morphology — more so than the 
merger history of all mergers. 

One might be tempted to conclude that the total major 
merger rate in the middle panel (major stellar mergers) 
is sufficiently steep to account for the change in mor- 
phological fraction with halo mass reported by W06 and 
W09; without bothering to account for the gas content 
of these mergers. Encouragingly, only ~ 30% of IO^^Mq 
halos have experienced a {m/M)staT > 0.3 merger since 
z = 2. It is important to note that even relatively minor, 
(m/A/)DM = 0.1 dark matter halo mergers with neg- 
ligible stellar content (m/M)star ^ 0.03 are capable of 
heating and thickening a galactic disk beyon d the prop- 
erties of the Milky Way (jPurcell et al.ll2009f ). (However, 
a galaxy need not contain a disk as thin as the Milky 
Way in order to be classified as late-type in either W06 or 
W09, so systems similar to those studied in lPurcell et al.l 
(|2009f l would still be labeled as "surviving" disks in our 
simple model.) Still, it is important to keep in mind that 
major dark matter mergers do not necessarily correspond 
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to major stellar mergers, esp ecially at halos less massiv e 
than the Milky Way (see e.g. lMalled[2008l : lStewart|[2009l ) . 
For example, consider a Mstar = lO^^M© galaxy experi- 
encing a stellar merger at 2; = 0, 1, 2 that is just below our 
definition of "major." with (m/M)star = 0.2. The corre- 
sponding dark matter halo ratios of these merger events 
will be ~ 0.4, 0.4, and 0.5, respectively. Even if we con- 
sider a more minor stellar merger, with (m/M)star = 0.1 
at these rcdshifts, such merger events still correspond to 
major dark matter mergers, with dark matter halo ratios 
of ^ 0.3, 0.3, 0.4, respectively. If gas content is ignored, a 
substantial fraction of these events will easily be capable 
of destroying disk morphologies altogether (not simply 
thickening and heating the existing disk). Still, these 
mergers would be classified as "minor" stellar mergers, 
with a stellar merger ratio < 0.3, even while the total 
dark matter mass ratios may approach 2:1. In addition, 
the opposite effect occurs at large galaxy masses. For ex- 
ample, a major stellar mergers into a Mgtar = IO^^Mq 
galaxy with a (m/M)star = 0.3 at z = only corresponds 
to a dark matter halo mass ratio of ~ 0.1, which is less 
likely to be morphologically destructive. This is why the 
total merger fraction in the middle panel (at high galaxy 
mass) exceeds that in the left panel: major stellar merg- 
ers only correspond to minor dark halo mergers at this 
mass regime. Thus, we conclude that the major stellar 
merger fractions in the middle panel likely present an 
uneven, and potentially incomplete picture of possible 
means of disk destruction via mergers. 

3.3. Gas Delivery Via Mergers 

In iStewart et all (|2009[ ). we discuss the observational 
implications of two well-known consequences of galaxy 
halo mergers: merger-induced starbursts and morpho- 
logical disturbance. A third potentially important con- 
sequence of mergers is the direct, cumulative deposition 
of cold baryons (gas and stars) onto galaxies. We are now 
concerned with the exact baryonic content of each merg- 
ing galaxy, whereas we have previously only focused on 
whether or not a given halo lies above or below an arbi- 
trary gas fraction threshold. As such, we impose an addi- 
tional constraint to our method for assigning gas (as out- 
lined in §2). Specifically, when estimating the gas content 
of halos at high redshift, blindly extrapolating equation 
[T]to arbitrarily small stellar masses sometimes results in 
more baryons in a given halo than the universal baryon 
fraction of matter in the Universe. This is an unphysical 
situation that arises from extrapolation far beyond the 
regime where Mgas(Mstar) is well-constrained. In order 
to avoid unphysically high gas content of low mass halos 
at high redshift, we present two models for assigning up- 
per limits to the gas content of low stellar mass galaxies. 
In the model A (left panel of figure [3|) , we set an upper 
limit on equation [T] such that Mgas + Mstav £ /b-A^DM, 
where fb = 0.17. In model B, we define the upper limit by 
the ratio of gas mass to halo mass: {fiim = Mgas/Muu 
at Mstar = 3 X 10^ Mq). For galaxies with stellar mass 
lower than this threshold, we then set Mgas = fumMDM- 
This model, which always assigns less (or equal) gas per 
galaxy than the first model, is used to construct the right 
panel of figure O 

In both panels of Figure [3l the solid black lines show 
the fraction of a central galaxy's current {z — 0) bary- 



onic mass (Mgai = Mgas + Mstar) acquired directly via 
major mergers as a function of halo mass /merged = 
Mi„crgcd/Mgai(z = 0). Specifically, Mmergo includes all of 
the baryonic mass in mergers obeying (m/M)DM > 0.3 
since z — 2. Focusing on the left panel (model A: 
Mgas + Mstar £ /&Mdm) the first clear result is that 
the merged baryonic fraction is significant: ~ 30 — 50% 
of the final galaxy mass is accreted directly in the form 
of major mergers. Given the effectiveness of dynami- 
cal friction in major mergers, we expect the majority 
of the accreted baryonic material in these events to be 
deposited in the central galaxy itself. In principle, this 
limits the fraction of a galaxy's baryons that can be ac- 
quired by direct hot halo cooling, cold flows, or minor 
mergers to 50 — 70%, depending on the halo mass of in- 
terest. Of course, (gaseous) baryons deposited via major 
mergers could in principle be blown out by energetic feed- 
back, but all in all, this result on major merger deposition 
(much like known results on cold flows) would seem to 
make the 'overcooling' problem in galaxy formation more 
difficult. While the fraction of baryons accreted directly 
as cold gas via major mergers is less substantial in the 
right panel (model B: Mgas < fumM-£,M), we still find 
significant accretion: 15 — 50% of the final galaxy mass. 

The dotted (blue) and dashed (red) fines in Figure [3] 
separate the total baryonic accretion fraction from major 
mergers into contributions from gas and stars, respec- 
tively (this is not a division between gas-rich and gas- 
poor mergers as before, but rather an integrated account- 
ing of all material regardless of the makeup of the merged 
progenitors). We see that in both models, the bary- 
onic accretion onto smaller halos (Mdm ^ 10^^ '^ Mq) 
is typically dominated by the gas content of the infalling 
galaxies, while the baryonic makeup of merged material 
into more massive systems is dominated by the infalling 
galaxies' stellar content. For Milky Way-size systems 
(Mdm = W^'^Mq), we find in model A (B) that typically 
~ 30% (20%) of a galaxy's baryonic content was accreted 
in the form of gas and stars directly via major mergers, 
with most of this accretion dominated by cold gas. In 
both models, more massive systems (Mdm = 10^^ Mq) 
typically accrete ~ 30%(10%) of their baryons as stars 
(gas) via major mergers, with no noticeable discrepancies 
between the two models. Though not shown, we find that 
most baryonic accretion from major mergers (70 — 80% 
of stars, 50 — 70% of gas) occurred at later times. 

How do our results change if we include more minor 
mergers in our accounting? If we count up all of the bary- 
onic acquisition in mergers larger than (?tt./M)dm > 0.1, 
we find that the mass fraction accreted as stars is boosted 
by a factor of of ~ 1.5 from the panels shown, while 
the accreted gas is amplified by a factor of ~ 1.7, both 
roughly independent of halo mass. We caution, how- 
ever, that the importance of minor mergers in deliver- 
ing baryons to central galaxies is significantly less clear 
than it is with major mergers. While the baryons asso- 
ciated with major mergers almost certainly become de- 
posited directly onto the central galaxy (see e.g. the 
simulations of Purcell et al. 2008b) the ultimate fate 
of the baryons in minor mergers will depend sensitively 
on the orbital properties of the secondary and on the 
potential presence of hot gas halo around the primary 
galaxy. Past work has demonstrated that the stellar ma- 
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Fig. 3. — The fraction of a 2 = galaxy's total baryonic mass that was accreted directly via major mergers or moderately sized mergers 
since 2: = 2, as a function of halo mass. Left: model A, in which we assign gas by Equation [T] but impose an upper limit due to the 
universal baryon fraction, such that Mgas + Mstar < fb^UM- Right: model B, in which we assign gas by Equation [T] but impose an upper 
limit such that Afgas < fiim^'^DM> where fujn at each redshift is the value of Mgaa/Mjjivi when Mgtar = 3 X 10* Mq at that redshift. In 
both panels, the dotted (blue) and dashed (red) lines show the accreted baryonic mass fraction from gas and stars, respectively, while the 
solid (black) line shows the total. 



terial in minor mergers will likely contribute to extended 
diffuse light co mponents like stellar halo s or intraclus- 
ter light fe.g.. iBullock fc JohnstoiJ l20()5l: iPiircell et all 
[2007HConrov et all I2007t IPurcell et al.ll2008D . but the 
destiny of accreted gas (which is the dominant compo- 
nent for galaxy halos) in these minor mergers is relatively 
unexplored. One possibility is that the gas in minor 
mergers i s quickly liberated v ia ram pressure stripping 
(see, e.g. iGrcevich et al.l [20081) and that it either evap- 
orates into the hot halo itself or eventually rains down 
onto the galaxy, possibly in the form of high-velocity 
clouds. These interesting possibilities are clearly beyond 
the scope of the present work but provide important av- 
enues for future investigation. 

The relation between gas mass and stellar mass below 
-^'^star = 3 X 10^ Mq remains relatively uncertain, making 
it impossible to know which panel of this figure is a more 
accurate representation of galaxy formation. In detail, 
we expect our two models to bracket reasonable expec- 
tations for the true relation between Mgtar and Mgas in 
this regime. Nevertheless, focusing on model A for the 
time being (left panel), we find it interesting that there 
appears to be a minimum in merged baryon fraction at 
the Mdm ^ IO^^Mq scale, which corresponds closely to 
the well-known mass scale of maximum galaxy forma- 
tion efficiency (^ in the galaxy luminosity function). 
Although it is unclear whether this minimum exists in 
reality, or is merely an artifice of the manner in which 
we have assigned gas, we speculate on a possible corre- 
spondence between galaxy formation efficiency and the 
merged baryon fraction. Specifically, a minimum in the 
merged baryon fraction naturally implies a maximum in 
the baryon fraction accreted via smooth gas accretion 
from the hot halo or from cold streams. We speculate 



that smooth gas accretion might allow for a higher ef- 
ficiency in star formation than accreting large clumps 
of gas via major mergers, because gas-rich mergers are 
likely to trigger massive starbursts that may blow sig- 
nificant gas content out of the central galaxy. We also 
speculate that accretion of stars via mergers may also 
be inefficient, since some fraction of a satellite galaxy's 
stars is likely to be distributed into the stellar halo before 
reaching the central galaxy within massive halos. If this 
is the case, that smooth gas accretion forms stars most 
efficiently, then it would be reasonable to expect a max- 
imum in the baryon fraction of smoothly accreted gas to 
correlate with the maximum galaxy formation efficiency. 
However, we emphasize that this possible correlation be- 
tween minimum merged baryonic fraction and maximum 
star formation efficiency is not a robust prediction of our 
model, but merely a speculation (for example, model B 
results in no such minimum in the merged baryonic frac- 
tion). 

3.4. Redshift Evolution 

While the cumulative fraction of halos that have ever 
experienced a gas-rich or gas-poor major merger (since 
z = 2) is the most pertinent question for morphologi- 
cal evolution and disk survivability (see §3.1), another 
point of interest is the redshift evolution of a more in- 
stantaneous measure of the merger rate of gas-poor and 
gas-rich mergers. Figure |4] shows the fraction of halos 
that have experienced at least one major merger with 
(m/M)DM > 0.3 in the past halo dynamical time, r ^. 

''As in lStewart et al.l I I2009I ) , in which we studied the evolution of 
the halo merger rate with redshift, we again adopt t{z) = R/V oc 
{A^{z) pu{z))~^/'^ , such that the halo dynamical time evolves with 



9 



„ 10.0 10.5 10.8 



11.1 



log Msior [Mg] 
10.0 10.5 10.8 1 1.1 



10.0 10.5 10.8 11.1 



o 



O 
C 

'(/) 



CP 



0) 



c 
o 



0.10 



t:o.oi 



z = 0.0 All Mergers 

T~1 .9 Gas Poor — — — - 
Gos Rich 

(m/M,)o„>0.3 




11.5 12.0 12.5 13.0 



11.5 12.0 12.5 13.0 
log [Mq] 



11.5 12.0 12.5 13.0 



Fig. 4. — Major merger fraction within the past dynamical time of the halo, t, as a function of halo mass. The solid (black) line shows the 
total merger fraction for dark matter halos (no baryons included). The dashed (red) line shows the merger fraction while only considering 
gas poor mergers {fg < 50%). The dotted (blue) line shows only gas-rich mergers (/g > 50%). The three panels show results for 2 = 0, 
z = 0.5, and 2 = 1, for which the halo dynamical time, r ~ 1.9, 1.3, 0.9 Gyr, respectively. We prim arily focus on the decomposition of these 
merger fractions into gas-rich and gas-poor mergers. We refer the reader to lStewart et all l|2009l ) for a detailed analysis of the dependence 
of dark matter of merger rates and merger fractions on redshift, halo mass, and mass ratio. 

As in Figure [2l the solid (black) line shows the total 
merger fraction for dark matter halos, while the dotted 
(blue) and dashed (red) lines show the major merger frac- 
tion for gas-poor {fg < 50%) and gas-rich {fg > 50%) 
mergers, respectively. Because the quantitative values of 
these merger fractions depend sensitively on the merger 
timescale in question at each redshift, we choose to fo- 
cus on two important qualitativ e resul t s from this figure. 
We refer the reader to [Stewart et all (|2009f ) for a more 
detailed discussion of the evolution of the halo merger 
rate (and merger fractions) with redshift, halo ma ss, and 
merger mass ratio (see also iFakhouri fc Mal[2009l for the 
effects of halo environment on the merger rate) . 

The first feature of note in this figure is the presence 
of a typical transition mass, {Mt)-DM, such that most of 
the recent major mergers into halos less massive than 
(Mt)DM are gas-rich, while most of the recent major 
mergers into halos more massive than {Mt)-DM are gas- 
poor. The existence of this transition mass is primarily 
due to the strong dependence of galaxy gas fractions on 
stellar mass (and thus, halo mass) at fixed redshift. In 
addition to creating this transition mass, the dependence 
of gas fraction on halo mass also results in a very limited 
mass range for which mixed mergers (where one galaxy 
is gas-poor and the other is gas-rich) constitute a sig- 
nificant portion of all mergers (at most ~ 50%). This 
effect is apparent in Figure [4] where the combined total 
of the gas-poor and gas-rich merger fractions fall signif- 
icantly short of the total. Not unexpectedly, this range 
of importance for mixed mergers is centered on (Mt)DM- 

The other important result from Figure [4] is that 
(Ait)DM is more massive at higher redshifts, with 

redshift, but is independent of halo mass. 



(Mt)DM ^ 10"■^ 10"■^ 1O12-«M0 at z = 0.0, 0.5, 1.0, re- 
spectively. This arises naturally from the strong increase 
in galaxy gas fractions to higher redshift. The corre- 
sponding galaxy stellar mass transitions at z = 0.0, 0.5, 
and 1.0, are (Mt)star 10^ ^ lO^O-^, 10" °Mo (upper 
horizontal axis in Figure 4). Of course, the precise value 
of Mt at each redshift will depend to some degree on 
our definitions of "gas-rich" and "gas-poor," but the ex- 
istence of this transition mass, and its qualitative evolu- 
tion with redshift should be robust to changes in these 
definitions. 

Consider recent major mergers into Milky- Way size 
IO^^A/q halos. At z = 0, mergers of this kind are very 
uncommon. Only ~ 5% of Milky- Way size halos should 
have experienced a major dark matter accretion event 
with (to/M)dm > 0.3 in the last r ~ 2 Gyr. However, 
when these major mergers do occur at z = they are very 
likely gas-poor (- 0.04/0.05 = 80% of the time). On the 
other hand major mergers are fairly common in IO^^Mq 
halos at z 1, with ~ 15% experiencing such a merger 
in the last r ~ 1 Gyr. Nevertheless, these higher redshift 
mergers are almost universally gas-rich. Under the pre- 
sumption that gas-rich mergers do not destroy disk mor- 
phologies, the evolution of the merger rate with redshift 
and in the associated gas-rich transition mass makes it 
increasingly likely that major mergers build disk galaxies 
at high redshift rather than destroy them (c.f. Robert- 
son et al. 2006a, Robertson & Bullock 2008). If we were 
to boldly extrapolate our trends to higher redshift, we 
would expect that nearly all major mergers into halos 
with A/dm < IO^^Mq should be gas-rich {fg > 50%) at 
z > 1. 

Encouraginglv. [Lin et al.l (|2008[ ) observe a similar red- 
shift evolution between gas-rich and gas-poor mergers 
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by studying the close-pair counts of galaxies from the 
DEEP2 Redshift Survey. While our definitions vary in 
detail from theirs (they divide galaxy pairs into wet and 
dry mergers based on galaxy colors, and bin their sample 
by total galaxy luminosity, while wc use galaxy gas frac- 
tions to define gas-poor versus gas-rich mergers, and bin 
by galaxy stellar mass) they also find that at fixed lu- 
minosity (stellar mass) the percentage of major mergers 
that are dry (gas-poor) should decrease with increasing 
redshift, while the percentage of major mergers that are 
wet (gas-rich) should increase. We reserve a more de- 
tailed comparison to their results for a future study, but 
we find the qualitative agreement encouraging. 

3.5. Comparison to Previous Work 

Recent studies of galaxy formation at high redshift 
using hydrodynamic simulations have stressed the im- 
portance of smooth accret ion of cold gas fro m filamen- 
tary streams. For example. iKeres et al.l ()2009[ ) compared 
the accretion rate of gas onto galaxies via cold flows 
and via mergers, and found that at z = 1 ~ 2, only 
about half of all gas accretion (onto galaxies correspond- 
ing to Mdm ^ IO^^'^Mq) is in the form of mergers, 
(where gas from mergers was defined as any gas that 
was added to galaxies i n dense baryonic clumps). Simi- 
larly, iDekel et all (|2009[ ) found that half of cold gas infall 
onto massive z — 2 galaxies is acquired via mergers with 
(m/A/)DM > 0.1, and the other acquired from cold flows. 
Indeed, even studies that focus primarily on the impor- 
tance of galaxy mergers also note that smooth gas ac- 
cretion is at least as dominant as galaxy mergers in the 
mass buildup of galaxies (e.g. iMaller et al.l 120061 ). Our 
results do not contradict these expectations. As demon- 
strated in Figure [31 wc expect that ^ 30 — 40% of a 
typical galaxy's baryons should have been accreted di- 
rectly via major mergers with (to/M)dm > 0.3. This 
leaves significant room for cold-flow gas to contribute to 
the baryonic assembly of small galaxies and for cooling 
to contribute to the buildup of larger galaxies. Though, 
as mentioned above, wc do expect that the percentage of 
merger-delivered baryons could rise to as much as ~ 60% 
if all of the baryonic material from (m/A/)DM > 0.1 
mergers is able to find its way into the central galaxy. 

iBrooks et al.l (|2009f ) used a high-resolution cosmolog- 
ical hydrodynamic simulation to study the gas accre- 
tion onto four disk galaxies within halos of masses 
-[^Qio.7-i2.7jy^^^ and also found that smooth accretion of 
gas (either shocked or un-shocked) dominates the mass 
buildup of their galaxies. When comparing smooth gas 
accretion to gas infall from mergers (using a generous 
definition of what qualifies as a merger) they found that 
~ 25% of the total gas infall into their Milky Way-size 
galaxy derives from mergers, with ~ 10% of the final 
stellar content at z = being accreted directly as stars 
from mergers. While their detailed results (and defini- 
tions) differ slightly from our own, the rough consistency 
between their simulation and our own semi-empirical ap- 
proach is quite encouraging. 

One point of caution associated with the discussion 
of cold flows is that these predictions arc based entirely 
on simulations that do not generally reproduce the ob- 
served baryonic mass function and stellar mass function 
of galaxies. It is possible that the cold flows are some- 
how restricted in the real universe in a way that solves 



the well-known over-cooling problem in galaxy forma- 
tion. Due to the inherent difficulties in detecting cold 
filaments of gas locally and at high redshift, there has 
yet to be an observational confirmation of a star forming 
galaxy fueled by the smooth accretion of cold gas along 
filaments, as seen in hydrodynamic simulations. In con- 
trast, our predictions for the accretion of stars and gas 
via major mergers is solidly normalized against observa- 
tions, and is arguably inevitable in the context of LCDM 
merger histories. Of course, baryonic material (especially 
gas) that is delivered via major mergers need not remain 
in the central galaxy indefinitely. Gas accreted either 
along cold flows or through major mergers may be sub- 
sequently expelled via supernovae or AGN feedback (e.g. 
Benson et al]|2003HDi Matteo et al.]|2005l : lSpringel et al.l 
2005; So merville et al.ll2008D . In this respect, Figure H 
represents an upper limit on the baryonic contribution 
from major mergers, with respect to the total amount of 
barons currently exist in the galaxy. 

4. CONCLUSION 

We have used dark matter halo merger trees from 
a large cosmology A^-body simulation together with 
observationally-normalized relationships between dark 
matter halo mass, galaxy stellar mass, and galaxy gas 
mass to explore the baryonic content of galaxy mergers 
back to redshift z — 2. Though our adopted associations 
between halo mass and the baryonic content of galaxies 
cannot be precisely correct, it is almost certainly accu- 
rate in its scalings with halo mass and redshift, and has 
the added advantage that it is independent of any uncer- 
tain galaxy formation physics. Indeed, any self-consistent 
galaxy formation model that is set within the LCDM 
framework would certainly need to reproduce our gross 
baryonic assignments in order to reproduce the observed 
universe. Our main results based on this methodology 
may be summarized as follows: 

1. The vast majority (^ 85%) of the major merg- 
ers experienced by Milky- Way size galaxies since 
z = 2 should have been gas-rich, and this frac- 
tion drops significantly towards higher mass sys- 
tems (see Figure 2). Remarkably, the fraction of 
galaxies with gas-poor major mergers matches well 
to the observed fraction of bulge-dominated galax- 
ies as a function of halo mass from A/dm ~ 10^^ to 

2. Though recent major mergers are expected to be 
rare for small galaxies in the local universe, the 
recent mergers that do occur should typically be 
gas poor. At higher redshift, recent mergers be- 
come more common and the probability that such 
a merger is gas-rich also increases (see Figure 4). 
One can define a transition dark matter halo mass 
Mt, below which most of the recent major merg- 
ers are gas-rich and above which they are gas 
poor, and this transition mass increases with red- 
shift: (M()dm ~ 10"'', 10"-^ 10^2.8^^ at z = 
0.0, 0.5, 1.0. As a result, the vast majority of recent 
major mergers into galaxy-size Mdm < IQ'^'^Mq 
dark matter halos are expected to be gas-rich at 
z<l. 



11 



3. A significant fraction (20 — 50%) of the baryonic 
mass in field galaxies at z = should have been 
deposited directly via major mergers since z = 2. 
For less massive galaxies, A/dm ~ 10^-^'^ Mq, the 
vast majority of the merger-acquired baryons are 
gaseous, while in more massive galaxies Mdm ~ 
IO^^Mq, major mergers bring in mostly stars (see 
Figure 3). For Milky Way-size systems, major 
mergers since z = 2 bring in 30% of the galaxy's 
z = baryonic mass, with most of this contribution 
in the form of gas. 

Many of these conclusions lend support to the con- 
j ecture of Robertson et al. (2006a) and iBrook et al.l 
(|2007ar i. who were the first to forcefully suggested a sce- 
nario where gas-rich mergers play an important role in 
building and stabilizing disk galaxies at high redshift. 
Though our conclusions are far from a sufficient test of 
this idea, we have demonstrated that gas-rich mergers 
should be common enough to make it viable for serious 
consideration. 

Among our most interesting results is the similarity 
between our predicted gas-poor merger fraction with 
halo mass and the observed early-type galaxy fraction 
with halo mass (Figure 2, left panel). Of course, even 
if gas-rich mergers do preserve disks, there are many 
openings for concern. For example, the current pre- 
sentation leaves little room for the production of bulge- 
dominated systems by means other than maj or mergers. 
In an extreme yet illustrative example, B ournaud et al.l 
(|2007l ) used a suite of focused simulations to show that 
bulge-dominated galaxies may be formed by successive 
minor mergers. Disk galaxies can also grow massive 
bulges by secular processes (typically bulges formed in 
this way show kinematically distinct properties from 
class ical bulges, and are referred to as "pseudobulges' 
e.g.lCourteau et aL|ll996 : iKorinendv fc Kennicuttl[200 



IKormendv fc Fisheil 12001 Hool 

An interesting possibility in this context of disk sur 
vival and secular evolution is that we have been too con 
servative in our classification of 'gas-rich'. Our fiducial 



division between gas-rich and gas poor at fg 



50% 



was motivated by the idealized simulations studied by 
Robertso n et al. (2006a) and Hopkins et al. (2008). 
However iGovernato et al.l (|2009f ) used a cosmologically 
self-consistent hydrodynamic simulation to demonstrate 
the creation of spiral galaxy at z = within a sys- 
tem that experienced a very major ((to/M)dm > 0.8) 
merger at z = 0.8. The two progenitor galaxies in this 



case were only moderately gas-rich {fg ~ 20%). Despite 
these relatively low gas fractions, the merger remnant 
was able to quickly reform a disk via the cooling of gas 
from the hot phase. If we use this result as motivation 
to focus on the more lenient {fg > 30%) definition of 
gas- rich in Figure [2l our gas-poor merger fractions drop 
to 10 — 20% smaller than the observed bulge-dominated 
fractions, leaving room for processes other than gas-poor 
major mergers to cause a significant portion of morpho- 
logical transformations. 

The general semi-empirical findings we have presented 
here may be regarded as accurate (not precise) predic- 
tions based on merger histories of LCDM halos and ob- 
served relations. As such, it is reassuring that our al- 
most unavoidable qualitative trends are consistent with a 
growing body of work that stresses the importance of gas- 
richness in preservi ng disk morphologies during mergers 



20051 Robertson et all l2006al: IBrook et all 


2007blla: 


Governato et al. 2007, 20091: iHookins et al. 


12009a: 


Robertson & BuUockl I2008D . Although we have focused 



primarily on issues of morphological transformation, disk 
survival, and baryonic accretion via mergers in this pa- 
per, we believe that in future work, the semi-empirical 
approach we have used here may provide a useful tool 
in exploring a vast array of galaxy properties and evolu- 
tionary mechanisms. 
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